March '93 



PUPT-1383 



LOOP EQUATION IN TURBULENCE 



A. A. Migdal 

Physics Department, Princeton University, 
Jadwin Hall, Princeton, NJ 08544-1000. 
E-mail: migdaWacm.princeton. edu 



Abstract 

The incompressible fluid dynamics is reformulated as dynamics of closed loops C 
in coordinate space. This formulation allows to derive explicit functional equation 
for the generating functional \&[C] in inertial range of spatial scales, which allows 
the scaling solutions. The requirement of finite energy dissipation rate leads then 
to the Kolmogorov index. We find an exact steady solution of the loop equation in 
inertial range of the loop sizes. The generating functional decreases as exp ^— 
where A = § c r A dr is the area inside the loop. The pdf for the velocity circulation 
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r is Lorentzian, with the width T oc A3 . 
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1 Introduction 

Incompressible fluid dynamics underlies the vast majority of natural phenomena. It is 
described by famous Navier-Stokes equation 

i) a = vdpv a - vpdpv a - d a p; d a v a = (1) 

which is nonlinear, and therefore hard to solve. This nonlinearity makes life more inter- 
esting, though, as it leads to turbulence. Solving this equation with appropriate initial 
and boundary conditions we expect to obtain the chaotic behavior of velocity field. 

The simplest boundary conditions correspond to infinite space with vanishing velocity 
at infinity We are looking for the translation invariant probability distribution for velocity 
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field, with infinite range of the wavelengths. In order to compensate for the energy 
dissipation, we add the usual random force to the Navier-Stokes equations, with the short 
wavelength support, corresponding to large scale energy pumping. 

One way to attempt to describe this probability distribution by the Hopf generat- 
ing functional (the angular bracket denote time averaging, or ensemble averaging over 
realizations of the random forces) 

Z[J] = ^exp (| d 3 rJ a (r)v a (r)^ (2) 

which is known to satisfy linear functional differential equation 

Z (3) 

similar to the Schrodinger equation for Quantum Field Theory, and equally hard to solve. 
Nobody managed to go beyond the Taylor expansion in source J , which corresponds to 
the obvious chain of equations for the equal time correlation functions of velocity field in 
various points in space. The same equations could be obtained directly from Navier-Stokes 
equations, so the Hopf equation looks useless. 

In this work we argue, that one could significantly simplify the Hopf functional without 
loosing information about correlation functions. This simplified functional depends upon 
the set of 3 periodic functions of one variable 

C:r a = C a (6); < 9 < 2tt (4) 

which set describes the closed loop in coordinate space. The correlation functions reduce 
to certain functional derivatives of our loop functional with respect to C{6) at vanishing 
loop C -> 0. 

The properties of the loop functional at large loop C also have physical significance. 
Like the Wilson loops in Gauge Theory, they describe the statistics of large scale structures 
of vorticity field, which is analogous to the gauge field strength. As we argue in this paper, 
the Kolmogorov scaling law corresponds to the loop functional decreasing as exp (—A%^j, 
where A is the tensor Area inside the loop. This area law emerges as a self-consistent 
solution of our loop equation in the inertial range of loops. By Fourier transformation 
of the loop functional we obtain the pdf for the velocity circulation, which turns out 
Lorentzian. 

In Appendix A we recover the expansion in inverse powers of viscosity by direct iter- 
ations of the loop equation. 

In Appendix B we study the matrix formulation of the Navier-Stokes equation, which 
may serve as a basis of the random matrix description of turbulence. 

In Appendix C we study the reduced dynamics, corresponding to the functional Fourier 
transform of the loop functional. We argue, that instead of 3D Navier-Stokes equations 
one can use the ID equations for the Fourier loop P a (9, t). 



Z = H 
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In Appendix D we discuss the relation between the initial data for velocity field and 
the P field, and we find particular realisation for these initial data in terms of the gaussian 
random variables. 

In Appendix E we introduce the generating functional for the scalar products P a (9)P a {9'). 
The advantage of this functional over the original ^[C] functional is the smoother con- 
tinuum limit. 

Finally, in Appendix F we discuss the possible numerical implementations of the re- 
duced loop dynamics. 

These four last Appendixes can be skipped at first reading. They might be needed for 
further development of this approach. 

2 The Loop Calculus 

We suggest to use in turbulence the following version of the Hopf functional 

= ^exp (i j> dC a {9)v a {C{9))^ (5) 

which we call the loop functional or the loop field. It is implied that all angular variable 9 
run from to 2n and that all the functions of this variable are 2ir periodic.^] The viscosity 
v was inserted in denominator in exponential, as the only parameter of proper dimension. 
As we shall see below, it plays the role, similar to the Planck's constant in Quantum 
mechanics, the turbulence corresponding to the WKB limit v — > 0. Q 

As for the imaginary unit i , there are two reasons to insert it in the exponential. 
First, it makes the motion compact: the phase factor goes around the unit circle, when 
the velocity field fluctuates. So, at large times one may expect the ergodicity, with well 
defined average functional bounded by 1 by absolute value. Second, with this factor of 
i , the irreversibility of the problem is manifest. The time reversal corresponds to the 
complex conjugation of ^, so that imaginary part of the asymptotic value of \P at t — ► oo 
measures the effects of dissipation. 

The loop orientation reversal C{9) — > C(2n — 9) also leads to the complex conjugation, 
so it is equivalent to the time reversal. This symmetry implies, that any correlator of 
odd/even number of velocities should be integrated odd/even number of times over the 
loop, and it must enter with an imaginary /real factor. Later, we shall use this property 
in the area law. 

1 This parametrization of the loop is a matter of convention, as the loop functional is parametric 
invariant. 

2 One could also insert any numerical parameter in exponential, but this factor could be eliminated by 
space- and/or time rescaling. 
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We shall often use the field theory notations for the loop integrals, 



* [C] = (exp ( - f dr a v, 



\v Jc 



(6) 



This loop integral can be reduced to the surface integral of vorticity field 



v. 




(7) 



by the Stokes theorem 




dr a v, 



L 



da^u^; dS = C 
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This is the well-known velocity circulation, which measures the net strength of the 
vortex lines, passing through the loop C. Would we fix initial loop C and let it move 
with the flow, the loop field would be conserved by the Euler equation, so that only the 
viscosity effects would be responsible for its time evolution. However, this is not what 
we are trying to do. We take the Euler rather than Lagrange dynamics, so that the loop 
is fixed in space, and hence $ is time dependent already in the Euler equations. The 
difference between Euler and Navier-Stokes equations is the time irreversibility, which 
leads to complex average ^> in Navier-Stokes dynamics. 

It is implied that this field \I/ [C] is invariant under translations of the loop C{9) — > 
C{6) + const. The asymptotic behavior at large time with proper random forcing reaches 
certain fixed point, governed by the translation- and scale invariant equations, which we 
derive in this paper. 

The general Hopf functional @) reduces for the loop field for the following imaginary 
singular source 



The \l/ functional involves connected correlation functions of the powers of circulation 
at equal times. 



This expansion goes in powers of effective Reynolds number, so it diverges in turbulent 
region. There, the opposite WKB approximation will be used. 

Let us come back to the general case of the arbitrary Reynolds number. What could 
be the use of such restricted Hopf functional? At first glance it seems that we lost most of 
information, described by the Hopf functional, as the general Hopf source J depends upon 
3 variables x,y, z whereas the loop C depends of only one parameter 9. Still, this infor- 
mation can be recovered by taking the loops of the singular shape, such as two infinitesimal 




(9) 




(10) 
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R- 



Fig. 1 



loops R\, i?2, connected by a couple of wires 
The loop field in this case reduces to 



vl/ [C] - (exp ( -£^uv( ri ) + — E>^(rs.) ) > (11) 



where 



E£, = jUr„r„ (12) 



is the tensor area inside the loop R. Taking functional derivatives with respect to the 
shape of R\ and R2 prior to shrinking them to points, we can bring down the product of 
vorticities at 7*1 and r 2 . Namely, the variations yield 

5YP V = f (drjr^ + r^d8r v ) = f (drjr^ - dr^5r v ) (13) 



IR JR 



where integration by parts was used in the second term. 

One may introduce th 
given point r at the loop. 



One may introduce the area derivative &a 5 ^ , which brings down the vorticity at the 
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v 



S 1 ^ [C] 



{uj^{ ri )u Xp {r 2 )) (14) 



The careful definition of these area derivatives are or paramount importance to us. 
The corresponding loop calculus was developed in[Q] in the context of the gauge theory. 
Here we rephrase and further refine the definitions and relations established in that work. 

The basic element of the loop calculus is what we suggest to call the spike derivative, 
namely the operator which adds the infinitesimal A shaped spike to the loop 

D a (9,e)= d4> l- 1 + ) 7777— r (15) 



5CM 
The finite spike operator 

A(r,9,e) = exp(r a D a (9,e)) (16) 
adds the spike of the height r. This is the straight line from C{9) to C(9+e)+r, followed by 
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C 



Fig. 2 



another straight line from C(9+e)+r to C(9+2e), 

Note, that the loop remains closed, and the slopes remain finite, only the second deriva- 
tives diverge. The continuity and closure of the loop eliminates the potential part of 
velocity; as we shall see below, this is necessary to obtain the loop equation. 

In the limit e — * these spikes are invisible, at least for the smooth vorticity field, as 
one can see from the Stokes theorem (the area inside the spike goes to zero as e) . However, 
taking certain derivatives prior to the limit e^Owe can obtain the finite contribution. 

Let us consider the operator 



By construction it inserts the smaller spike on top of a bigger one, in such a way, that a 



polygon appears 

Taking the derivatives with respect to the vertices of this polygon r, r' , setting r = r' = 
and antisymmetrising, we find the tensor operator 




(17) 




r 



Fig. 3 




(18) 



which brings down the vorticity, when applied to the loop field 



Q a/3 (9,e)^[C] 



u} aP (c(e))*[C\ 



(19) 
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The quick way to check these formulas is to use formal functional derivatives 

Taking one more functional derivative derivative we find the term with vorticity times 
first derivative of the 5 function, coming from the variation of C'(6) 

5 2 9\C] 59 \C] S 2 9\C] 



8C a {9)8C p {6') y J 5a af3 (C(9)) ^' AV ' 5a ai (C(8)) 5a px (C(9>)) 

This term is the only one, which survives the limit e — > in our relation (|T^) . 

So, the area derivative can be defined from the antisymmetric tensor part of the second 
functional derivative as the coefficient in front of 5'(9 — 9') . Still, it has all the properties 
of the first functional derivative, as it can also be defined from the above first variation. 
The advantage of dealing with spikes is the control over the limit e — > , which might be 
quite singular in applications. 

So far we managed to insert the vorticity at the loop C by variations of the loop field. 
Later we shall need the vorticity off the loop, in arbitrary point in space. This can be 
achieved by the following combination of the spike operators 

A{r,e,e)U{r 1 ,r 2 ,e + e,8); 5< e (22) 

This operator inserts the II shaped little loop at the top of the bigger spike, in other 
words, this little loop is translated by a distance r by the big spike. 

Taking derivatives, we find the operator of finite translation of the vorticity 

A(r,9,e)Q a p(9 + e,5) (23) 
and the corresponding infinitesimal translation operator 

D^(9,e)n al3 (9 + e,5) (24) 
which inserts d^Uap (C(9)) when applied to the loop field. 

Coming back to the correlation function, we are going now to construct the operator, 
which would insert two vorticities separated by a distance. Let us note that the global A 
spike 

AM,,)=exp( r „f (25) 
when applied to a shrunk loop C(<f>) = does nothing but the backtracking from to r 

< > r 

Fig. 4 

This means that the operator 

n Ql3 (0,5)n Xp (7T,6)A(r,0,7T) (26) 
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when applied to the loop field for a shrunk loop yields the vorticity correlation function 

n aP (0, 5)n Xp (ir, S)A (r, 0, tt) [0] = (u; af3 (0)uj Xp (r)) (27) 

The higher correlation functions of vorticities could be constructed in a similar fashion, 
using the spike operators. As for the velocity, one should solve the Poisson equation 



with the proper boundary conditions , say, v = at infinity. Formally, 

v a (r) = ^dpu Pa {r) 

This suggests the following formal definition of the velocity operator 

1 



(28) 



(29) 



V a (e,e,6) 



D p (9,e)n Pa (9,S); 5 « e 



V a (9,e,5MC] S ^°v a (C(9))nC] 



(30) 
(31) 



Another version of this formula is the following integral 

V a (9, e,S)=J d 3 p^^A (p, 9, e) Q a/3 (9 + e, 6) (32) 
where the A operator shifts the £1 by a distance p off the original loop at the point r = 





Fig. 5 



C{9+e) 



3 Loop Equation 

Let us now derive exact equation for the loop functional. Taking the time derivative of 
the original definition, and using the Navier-Stokes equation we get in front of exponential 



jf drj- (vdpv a - vpdpVa - d a p) 



(33) 
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The term with the pressure gradient yields zero after integration over the closed loop, and 
the velocity gradients in the first two terms could be expressed in terms of vorticity up 
to irrelevant gradient terms, so that we find 

j> dr a - {vdpupa - vpwpa) (34) 

Replacing the vorticity and velocity by the operators discussed in the previous Section 
we find the following loop equation (in explicit notations) 

-t*[C]= (35) 

j>dC a {9) (Dp(e,e)np a {9,e) + ^J ^p^jgA (p, 9, e) Q^(9 + e,5)n Pa {9, 5)^ *[C] 

The more compact form of this equation, using the notations of , reads 

i uif[C] = H C V (36) 



lc \ bop a {r) J An\r - r'\ 6 da l 3 a {r)bap~ j { r ) J 
Now we observe that viscosity v appears in front of time and spatial derivatives, like the 
Planck constant H in Quantum mechanics. Our loop hamiltonian Tic is not hermitean, 
due to dissipation. It contains the second loop derivatives, so it represents a (nonlocal!) 
kinetic term in loop space. 

So far, we considered so called decaying turbulence, without external energy source. 
The energy 

E = j d\\ vl (37) 

would eventually all dissipate, so that the fluid would stop. In this case the loop wave 
function \l/ would asymptotically approach 1 

*[C] ^ 1 (38) 

In order to reach the steady state, we add to the right side of the Navier-Stokes 
equation the usual gaussian random forces f a (r,t) with the space dependent correlation 
function 

(f a (r, t)f p {r', 0) = M(* - Wr ~ r') (39) 
concentrated at at small wavelengths, i.e. slowly varying with r — r'. 
Using the identity 

(/ a (r, *)*[«(.)]> = J d r ' F ( r - r ')j^jf ( 4 °) 

which is valid for arbitrary functional $ we find the following imaginary potential term 
in the loop hamiltonian 

SH C = % U[C] = - I dr a I dr' a F(r - r') (41) 



Note, that orientation reversal together with complex conjugation changes the sign of 
the loop hamiltonian, as it should. The potential part involves two loop integrations times 
imaginary constant. The first term in the kinetic part has one loop integration, one loop 
derivative times imaginary constant. The second kinetic term has one loop integration, 
two loop derivatives and real constant. The left side of the loop equation has no loop 
integrations, no loop derivatives, but has a factor of i . 

The relation between the potential and kinetic parts of the loop hamiltonian depends 
of viscosity, or, better to say, it depends upon the Reynolds number, which is the ratio 
of the typical circulation to viscosity. In the viscous limit, when the Reynolds number 
is small, the loop wave function is close to 1. The perturbation expansion in - goes in 
powers of the potential, in the same way, as in Quantum mechanics. The second (nonlocal) 
term in kinetic part of the hamiltonian also serves as a small perturbation (it corresponds 
to nonlinear term in the Navier-Stokes equation). The first term of this perturbation 
expansion is just 



with F(k) being the Fourier transform of F(r). This term is real, as it corresponds to the 
two- velocity correlation. The next term comes from the triple correlation of velocity, and 
this term is purely imaginary, so that the dissipation shows up. 

This expansion can be derived by direct iterations in the loop space as in ||, inverting 
the operator in the local part of the kinetic term in the hamiltonian. This expansion 
is discussed in Appendix A. The results agree with the straightforward iterations of the 
Navier-Stokes equations in powers of the random force, starting from zero velocity. 

So, we have the familiar situation, like in QCD, where the perturbation theory breaks 
because of the infrared divergencies. For arbitrarily small force, in a large system, the 
region of small k would yield large contribution to the terms of the perturbation expansion. 
Therefore, one should take the opposite WKB limit v — ► 0. 

In this limit, the wave function should behave as the usual WKB wave function, i.e. 
as an exponential 



The effective loop Action S[C] satisfies the loop space Hamilton- Jacobi equation 



The imaginary part of S[C] comes from imaginary potential U[C], which distinguishes 
our theory from the reversible Quantum mechanics. The sign of 35 must be positive 
definite, since \^\ < 1. As for the real part of 5[C], it changes the sign under the loop 
orientation reversal C(8) — > C(2n — 9). 




(42) 




(43) 




(44) 
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At finite viscosity there would be an additional term 

- v I dr a d,p^- -tvldrj rfV /;" r ; 3 fff" (45) 
Jc o<jp a {r) Jc J 47r|r — r'| J oa/3 a (r)dapy[r) 

on the right of (fHj). As for the term 

. dr a i {d p S[C\) (46) 
c oap a {r) 

which formally arises in the loop equation, this term vanishes, since dpS[C] = 0. This 
operator inserts backtracking at some point at the loop without first applying the loop 
derivative at this point. As it was discussed in the previous Section, such backtracking 
does not change the loop functional. This issue was discussed at length in 0, where the 
Leibnitz rule for the operator d a -^— was established 



iff = f'(9lC])a a ^L (47) 

In other words, this operator acts as a first order derivative on the loop functional with 
finite area derivative (so called Stokes type functional). Then, the above term does not 
appear. 

The Action functional S[C] describes the distribution of the large scale vorticity struc- 
tures, and hence it should not depend of viscosity. In terms of the above connected 
correlation functions of the circulation this corresponds to the limit, when the effective 
Reynolds number goes to infinity, but the sum of the divergent series tends to the 
finite limit. According to the standard picture of turbulence, the large scale vorticity 
structures depend upon the energy pumping, rather than the energy dissipation. 

This , of course implies, that both time t and the loop sizef] |C| should be greater then 
the viscous scales 

t > t = |C| > r = i/ff-i (48) 

where £ is the energy dissipation rate. 

It is defined from the energy balance equation 

° = &t (l^) = V ( Va ° 2va ) + ( faVa ^ (49) 
which can be transformed to 

^(u^) = 3F(0) (50) 

The left side represents the energy, dissipated at small scale due to viscosity, and the 
right side - the energy pumped in from the large scales due to the random forces. Their 
common value is S. 



! As a measure of the loop size one may take the square root of the minimal area inside the loop. 
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We see, that constant F(r — r'), i.e., F(k) oc 6(h) is sufficient to provide the necessary- 
energy pumping. However, such forcing does not produce vorticity, which we readily see 
in our equation. The contribution from this constant part to the potential in our loop 
equation drops out (this is a closed loop integral of total derivative). This is important, 
because this term would have the wrong order of magnitude in the turbulent limit - it 
would grow as the Reynolds number. 

Dropping this term, we arrive at remarkably simple and universal functional equation 



S[C] 



c 



dr n Id r 



5S 



5S 



An\r — r'\ 3 8ap a (r) 5a, 



(51) 



The stationary solution of this equation describes the steady distribution of the circulation 
in the strong turbulence. Note, that the stationary solutions come in pairs ±S. The sign 
should be chosen so, that 3$ > 0, to provide the inequality \^/\ < 1. 



4 Scaling law 

The 'Hamilton- Jacobi' equation without the potential term flBT] ) allows the family of the 
scaling solutions 



s[c] = 



c 



with arbitrary index k. The scaling function satisfies the equation 



(2«-1)0[C]-k/ dT a P^rr p = 

Jc oop a {r) ./<■ 



(52) 



(53) 



4n\r — r'| 3 8a/3 a (r) 5ap y (r') 

The left side here was computed, using the chain rule differentiation of functional. 

Asymptotically, at large time, we expect the fixed point, which is the homogeneous 
functional 



C 



\c\ 



S^C] = \C\ 2 ~f 
zeroing the right side of our 'kinetic' functional equation 

= 



/dV r',- r, 6SM SSM 



c 



47r|r — r'| 3 5ap a (r) 5aa^(r 



Tpy{ 



(54) 



(55) 



The Kolmogorov scaling would correspond to 

3 



K 



in which case one can express the S functional in terms of £ 

C 



S[C] = St 1 



(56) 



(57) 
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One can easily rephrase the Kolmogorov arguments in the loop space. The relation 
between the energy dissipation rate and the velocity correlator reads 

S = (v a {r o )v {O)d v a {O)) (58) 

where the point splitting at the viscous scale r is introduced. Such splitting is necessary 
to avoid the viscosity effects; without the splitting the average would formally reduce to 
the total derivative and vanish. 

Instead of the point splitting one may introduce the finite loop of the viscous scale 
\C\ ~ r , and compute this correlator in presence of such loop. This reduces to the WKB 
estimates 

*«*(r) - t^t; v a (r) = f dV /J'"] ; awCO (59) 
0Caj3{rj J 47r|r — r | d 

Using the generic scaling law for S we find 

_1 1-1 2-1 

u;~r K ;f~r K ; £ ~ r K (60) 

We see, that the energy dissipation rate would stay finite in the limit of the vanishing 
viscous scale only for the Kolmogorov value of the index. This argument looks rather 
cheap, but I think it is basically correct. The constant value of the energy dissipation 
rate in the limit of vanishing viscosity arises as the quantum anomaly in the field theory, 
through the finite limit of the point splitting in the correspondent energy current .0 

There is another version of this argument, which I like better. The dynamics of Euler 
fluid in infinite system would not exist, for the non- Kolmogorov scaling. The extra powers 
of loop size would have to enter with the size L of the whole system, like (^jpj ■ So, in 
the regime with finite energy pumping rate S the infinite Euler system can exist only for 
the Kolmogorov index. This must be the essence of the original Kolmogorov reasoning 

i- 

The problem is that nobody proved that such limit exists, though. Within the usual 
framework, based on the velocity correlation functions, one has to prove, that the infrared 
divergencies, caused by the sweep, all cancel for the observables. Within our framework 
these problems disappear, as we shall see later. 

As for the correlation functions in inertial range, unfortunately those cannot be com- 
puted in the WKB approximation, since they involve the contour shrinking to a double 
line, with vanishing area inside. Still, most of the physics can be understood in loop terms, 
without these correlation functions. The large scale behavior of the loop functional reflects 
the statistics of the large vorticity structures, encircled by the loop. 



4 I am grateful to A. Polyakov and E. Siggia for inspiring comments on this subject. 
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5 Area law 



The Wilson loop in QCD decreases as exponential of the minimal area, encircled by the 
loop, leading to the quark confinement. What is the similar asymptotic law in turbulence? 
The physical mechanisms leading to the area law in QCD are absent here. Moreover, there 
is no guarantee, that \I/[C] always decreases with the size of the loop. 

This makes it possible to look for the simple Anzatz, which was not acceptable in 
QCD, namely 

S[C] = s (E£) (61) 

where 

E° = I r^dr u (62) 



is the tensor area encircled by the loop C . The difference between this area and the 
scalar area is the positivity property. The scalar area vanishes only for the loop which 
can be contracted to a point by removal of all the backtracking. As for the tensor area, 
it vanishes, for example, for the 8 shaped loop, with opposite orientation of petals. 

Thus, there are some large contours with vanishing tensor area, for which there would 
be no decrease of the \& functional. In QCD the Wilson loops must always decrease at 
large distances, due to the finite mass gap. Here, the large scale correlations are known 
to exist, and play the central role in the turbulent flow. So, I see no reasons to reject the 
tensor area Anzatz. 

This Anzatz in QCD not only was unphysical, it failed to reproduce the correct short- 
distance singularities in the loop equation. In turbulence, there are no such singularities. 
Instead, there are the large-distance singularities, which all should cancel in the loop 
equation. 

It turns out, that for this Anzatz the (turbulent limit of the) loop equation is satisfied 
automatically, without any further restrictions. Let us verify this important property. 
The first area derivative yields 

n , , 5S ds , 

< w = s^vr 2 w, (63) 

The factor of 2 comes from the second term in the variation 

7 — = 8 ail 8p v - 5 av 5p^ (64) 
dcv(r) 

Note, that the right side does not depend on r. Moreover, you can shift r aside from the 
base loop C, with proper wires inserted. The area derivative would not change, as the 
contribution of wires drops. 

This implies, that the corresponding vorticity oo^ u (r) is space independent, it only 
depends upon the loop itself. The velocity can be reconstructed from vorticity up to 
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irrelevant potential terms 

v C p{r) = \r a u% (65) 
This can be formally obtained from the above integral representation 

as a residue from the infinite sphere R — \r'\ —>■ oo. One may insert the regularizing factor 
\r'\~ e in to, compute the convolution integral in Fourier space and check that in the limit 
e — > + the above linear term arises. So, one can use the above form of the loop equation, 
with the analytic regularization prescription. 

Now, the v u term in the loop equation reads 

£ dr 7 v c p (r) u% f oc u% to c Pl (67) 

This tensor trace vanishes, because the first tensor is antisymmetric, and the product of 
the last two antisymmetric tensors is symmetric with respect to cry. 

So, the positive and negative terms cancel each other in our loop equation, like the 
"income" and "outcome" terms in the usual kinetic equation. We see, that there is an 
equilibrium in our loop space kinetics. 

From the point of view of the notorious infrared divergencies in turbulence, the above 
calculation explicitly demonstrates how they cancel. By naive dimensional counting these 
terms were linearly divergent. The space isotropy lowered this to logarithmic divergency 
in (|66|) , which reduced to finite terms at closer inspection. Then, the explicit form of 
these terms was such, that they all cancelled. 

This cancellation originates from the angular momentum conservation in fluid mechan- 
ics. The large loop C creates the macroscopic eddy with constant vorticity u^g and linear 
velocity v c (r) cx r. This is a well known static solution of the Navier-Stokes equation. 
The eddy is conserved due to the angular momentum conservation. The only nontrivial 
thing is the functional dependence of the eddy vorticity upon the shape and size of the 
loop C. This is a function of the tensor area E£,, rather than a general functional of the 
loop. 

Combining this Anzatz with the space isotropy and the Kolmogorov scaling law, we 
arrive at the turbulent area law 

^ocexpf-B (^(Sy 2 )") (68) 

The universal constant B here must be real, in virtue of the loop orientation symmetry. 
When the orientation is reversed C(9) —>■ C(2n — 9), the loop integral changes sign, but 
its square, which enters here, stays invariant. Therefore, the constant in front must be 
real. The time reversal tells the same, since both viscosity v and the energy dissipation 
rate £ are time-odd. Therefore, the ratio ^ is time-even, hence it must enter \&[C] with 
the real coefficient. Clearly, this coefficient B must be positive, since |^[C]| < 1. 



15 



6 Discussion 



So, we found an exact solution of the loop equation in the turbulent limit. It remains to 
be seen, whether this is the most general solution, and is it realized in turbulent flows. 
Meanwhile, let us discuss its general properties, and its implications to the large scale 
vorticity distribution. 

First of all, let us address the issue of the uniqueness of this solution. Let us take the 
following Anzatz 

S[C\ = fUdr a £dr f a W(r-r f j\ (69) 

When substituted into the static loop equation (with the area derivatives computed in 
Appendix A), it yields the following equation for the correlation function W(r) 

= £ dr a £ dr'p £ dr'^U^r, r', r") (70) 
U a/ 3j(r, r', r") = W(r - r')V^W{r' - r") + permutations 

The derivative /' of the unknown function drops from the static equation. 

This equation should hold for arbitrary loop C. Using the Taylor expansion for the 
Stokes type functional ||, we can argue, that the coefficient function U must vanish up 
to the total derivatives. An equivalent statement is that the third area derivative of this 
functional must vanish. Using the loop calculus (see Appendix A) we find the following 
equation 

= V^'%V''^ v nU M ,{r, r', r") (71) 

which should hold for arbitrary r, r', r" . This leads to the overcomplete system of equations 
for W(r) in general case. However, for the special case W(r) = r 2 which corresponds to 
the square of the tensor area 

Kp = -\f c dr a j c dr' a {r-r'Y (72) 

the system is satisfied as a consequence of certain symmetry. In this case we find in the 
loop equation 

2 * dr a f dr'p(r - r') 2 * dr" a (r'p - r'p) oc j> dr a j> dr'^ir - r'f (73) 

The last integral is symmetric with respect to permutations of a, (3, whereas the first 
factor E?g is antisymmetric, hence the sum over af3 yields zero, as we already saw above. 

This solution can also be used to get the probability distribution for the velocity 
circulation. For this purpose the extra factor 7 should be inserted in the definition of the 
loop average 

*[Cl = (exp(^£dr a v a ^ (74) 
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and the Fourier transformation should be performed 



Pc(T) = -L / + °°rf 7 exp (--^f) *[C] (75) 



2ttz/ J- 

Now, the 7 dependence of ^f[C] can be found from the dimensional analysis. It enters 
only in combination with viscosity u, and the other source of viscosity, the vd 2 v term, 
was dropped in the Navier-Stokes equation. So, the asymptotic solution (|68| ) should get 
the extra factor of \^\ in the exponential. 

m oc exp ( -B | 7 | (^(^) 2 )"j (76) 
Another way to get this factor is to observe, that one could scale 7 away from original 



definition, by rescaling 8 — > 8\"f\ , since the velocity circulation scales as 8 
The Fourier integral yields the Lorentz distribution 



fc(r) = -^— T |f = B(f(ES,) ) (77) 



7T 

Let us answer an obvious question. The Lorentz distribution is symmetric. How does it 
agree with the known asymmetry of velocity correlations, in particular, the Kolmogorov 
triple correlation? The answer is that the Kolmogorov correlation does not imply the 
asymmetry of vorticity correlations. 

Taking the tensor version of the | law in arbitrary dimension d 

8 ( 2 \ 

and differentiating, we find that 

M0) V/3 (0V 7A (r)> = (79) 

So, the odd vorticity correlations could, in fact, be absent, in spite of the asymmetry 
of the velocity distribution. Besides, the area law does not apply to small loops which 
are involved in the definition of the vorticity correlations in terms the loop functional. 
Moreover, we do not see any reason to expect these correlations to be scale invariant. 
From our point of view these are not the asymptotic quantities, so the viscous effect 
could be important. 

With velocity correlations it is even worse than that. The infrared divergencies may 
be important as well, so that the factors like ( ^ could appear in the higher moments of 
velocity distribution. The observed violations of the Kolmogorov scaling in these moments 
could be attributed to these infrared divergencies. However, according to our theory, the 
velocity circulation agrees with the Kolmogorov scaling and has smooth Lorentz distribu- 
tion in the infinite system. 
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It would be extremely interesting to measure the velocity circulation for large loops 
in real or numerical experiments, and test these predictions. Maybe the long tails of the 
Lorentz distribution reflect the notorious intermittency? Note, that all the even moments 
are infinite for the Lorentz distribution, which is quite unusual, but does not contradict 
any known physical requirements. 

In real world this would probably mean that the n-th moments of the circulation are 
cut off by the finite time and finite size effects, i.e. they grow as min j(£t 2 ) n , (SL 4 ) 1 j. 
Instead of measurements of these moments, the shape of the Lorentz distribution and the 
area dependence of its width f could be tested. 

It was assumed in above arguments, that the loop C consist of only one connected 
part. Let us now consider the more general situation, with arbitrary number n of loops 
Ci, . . . C n . The corresponding Anzatz would be 

S n [C 1 ,...C n ] = s n (z\...Z n ) (80) 

where are tensor areas. 

This function should obey the same WKB loop equations in each variable. Introducing 
the loop vorticities 



which are constant on each loop, we have to solve the following problem. What are the 
values of io k such that the single velocity field v a (r) could produce them? 



We do not see any other solutions, but the trivial one, with all equal uA, and linear 
velocity, as before. This would correspond to 

s n (E\...E n ) = Sl (E); Z„„ = J2Kv= <f r >„ ( 82 ) 

k=i J<sCk 

The loop equation would be satisfied like before, with C = l±lC fe . This corresponds to the 
additivity of loops 

S n [C 1 ,...C n ] = S 1 [\HC k ] (83) 

Note, that such additivity is the opposite to the statistical independence, which would 
imply that 

S n [C 1 ,...C n ]='£S 1 [C k ] (84) 

The additivity could also be understood as a statement, that any set of n loops is equiva- 
lent to a single loop for the abelian Stokes functional. Just connect these loops by wires, 
and note that the contribution of wires cancels. So, if the area law holds for arbitrary 
single loop, than it must be additive. 

This assumption may not be true, though, as it often happens in the WKB approxi- 
mation. There is no single asymptotic formula, but rather collections of different WKB 
regions, with quantum regions in between. In our case, this corresponds to the following 
situation. 
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Take the large circular loop, for which the WKB approximation holds, and try to split 
it into two large circles. You will have to twist the loop like the infinity symbol oo, in 
which case it intersects itself. At this point, the WKB approximation might break, as 
the short distance velocity correlation might be important near the self- intersection point. 
This may explain the paradox of the vanishing tensor area for the oo shaped loop. From 
the point of view of our area law such loop is not large at all. 

This cancellation of large circulations is puzzling. Apparently, the further analytical 
and numerical study of the loop equation is required. One could try some variational 
approach, by approximating the loop by a polygon, which reduces the \l/ functional to a 
function of the vertices of the polygon. There are some other options, such as truncation 
of Fourier expansion of the loop function C{6). I do not expect any easy success here, 
because of singularities of the loop derivatives involved. Still, the beauty of the loop 
dynamics and its apparent reduction of the dimension of the turbulence problem, raises 
some hopes of the analytical advances. 

The main issue, in my opinion, is the Kolmogorov scaling. Usually, the indexes in the 
scaling solutions of the nonlinear integral equations of QFT (so called bootstrap equa- 
tions), are determined from selfconsistency of the equations, rather than from some extra 
requirements. However, the Kolmogorov law for triple correlation function was derived 
from the same Navier-Stokes equation plus the scaling assumptions about velocity corre- 
lation functions, which is equivalent to the assumptions we made in our scaling Anzatz 
for S[C]. So, it is possible, that the dynamical value of k, found from selfconsistency, 
would coincide with the Kolmogorov value. 

However, we have to keep in mind the possibility of the more general phenomenon. 
Namely, there might be the spectrum of solutions for the critical index k, corresponding 
to various fixed points of above Hamilton- Jacobi equation in the loop space. The steady 
state may not exist, if these fixed points are all unstable. In this case, S[C] would go 
from one fixed point to another. 

This would be the next level of complexity, as compared to the strange attractors, 
found in dynamical systems with small number of degrees of freedom. Not just the 
trajectory in phase space, but the whole probability distribution functional would evolve 
with with time. When averaged over given time interval T, it would reach steady state 
only for the scales less than T K . The marginal scales distribution would slowly drift, and 
larger scales would not be in equilibrium at all. In this case the infinite system would 
never reach the steady state. 

This situation is not described by above area laws, but rather requires more general 
non-steady solutions of the loop equations. In Appendices C, D, E, F we develop the 
general framework for studying such solutions. We decided to place these parts of our 
work in Appendix because this formalism is too heavy. However, the mathematically 
oriented reader might find it useful. 
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A Loop Expansion 



Let us outline the method of direct iterations of the loop equation. The full description 
of the method can be found in The basic idea is to use the following representation 
of the loop functional 

^] = l + E;((^-{<"} r W^ an ( ri ,...r n ) (85) 
„_ n Uc Jc ) cyclic 



cyclic 

This representation is valid for every translation invariant functional with finite area 
derivatives (so called Stokes type functional). The coefficient functions W can be related 
to these area derivatives. The normalization ^[0] = 1 for the shrunk loop is implied. 

In general case the integration points r 1; ...r n in (|85|) are cyclicly ordered around 
the loop C. The coefficient functions can be assumed cyclicly symmetric without loss 
of generality. However, in case of fluid dynamics, we are dealing with so called abelian 
Stokes functional. These functionals are characterized by completely symmetric coefficient 
functions, in which case the ordering of points can be removed, at expense of the extra 
symmetry factor in denominator 

= 1 + £ ^ I dr? ..J dr a n "W2^ an (n, . . . r n ) (86) 

n=2 C C 

The incompressibility conditions 

d ak W^ a J ri ,...r n ) = (87) 

does not impose any further restrictions, because of the gauge invariance of the loop 
functionals. This invariance (nothing to do with the symmetry of dynamical equations!) 
follows from the fact, that the closed loop integral of any total derivative vanishes. So, the 
coefficient functions are defined modulo such derivative terms. In effect this means, that 
one may relax the incompressibility constraints fl37|), without changing the loop functional. 

To avoid confusion, let us note, that the physical incompressibility constrains are not 
neglected. They are, in fact, present in the loop equation, where we used the integral 
representation for the velocity in terms of vorticity. Still, the longitudinal parts of XV 
drop in the loop integrals. 

The loop calculus for the abelian Stokes functional is especially simple. The area 
derivative corresponds to removal of one loop integration, and differentiation of the cor- 
responding coefficient function 

SrU = E JJ £ *? -£ «;...«. ru-r „) (88) 



where 



V° = d^S ua - cU MQ (89) 
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In the nonabelian case, there would also be the contact terms, with W at coinciding 
points, coming from the cyclic ordering 0. In abelian case these terms are absent, since 
W is completely symmetric. 

As a next step, let us compute the local kinetic term 

US![C\ = I dr^P^r (90) 
Jc oa^ v {r) 

Using above formula for the loop derivative, we find 

= E A / ^ I dr ^ ■■■I dr n nd2w ^l.a n (r,n,... r n ) (91) 
n= 1 ^* J C JC J c 

The net result is the second derivative of W with respect to one variable. Note, that the 
second term in V* dropped, as the total derivative in the closed loop integral. 

As for the nonlocal kinetic term, it involves the second area derivative off the loop, 
at the point r', integrated over r' with the corresponding Green's function. Each area 
derivative involves the same operator V, acting on the coefficient function. Again, the 
abelian Stokes functional simplifies the general framework of the loop calculus. The 
contribution of the wires cancels here, and the ordering does not matter, so that 

5 2 V[C] ~ 1 



5<y^{r)5(r^ vl {r') grille " Jc 



E ry f dr ? ■■■{ d < n K^ W Zl^ n r', r l5 . . . r n ) (92) 

„ n 7v* JC JC 



Lb ( w 1 



n+1 



Using these relations, we can write the steady state loop equation as follows 
Here the light dotted lines symbolize the arguments aifc, of W, the big circle denotes the 
loop C, the tiny circles stand for the loop derivatives, and the pair of lines with the arrow 
denote the Green's function. The sum over the tensor indexes and the loop integrations 
over rfc are implied. 

The first term is the local kinetic term, the second one is the nonlocal kinetic term, 
and the right side is the potential term in the loop equation. The heavy dotted line in 
this term stands for the correlation function F of the random forces. Note that this term 
is an abelian Stokes functional as well. 
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The iterations go in the potential term, starting with \1/[C] = 1. In the next approx- 
imation, only the two loop correction W 2 l0l2 (r\,r2) is present. Comparing the terms, we 
note, that nonlocal kinetic term reduces to the total derivatives due to the space symme- 
try (in the usual terms it would be (vco) at coinciding arguments), so we are left with the 
local one. 



This yields the equation 



v 3 d 2 W^(r-r') = F(r-r>)5 a p 



modulo derivative terms. The solution is trivial in Fourier space 

W^(r - r ') = -J^ exp (, k(r - r')) 6 a / {k) 



(2nY 



Note, that we did not use the transverse tensor 



v 3 k 2 



(93) 



(94) 



(95) 



Though such tensor is present in the physical velocity correlation, here we may use 5 a p 
instead, as the longitudinal terms drop in the loop integral. This is analogous to the 
Feynman gauge in QED. The correct correlator corresponds to the Landau gauge. 

The potential term generates the four point correlation FW 2 . which agrees with the 
disconnected term in the W 4 on the left side 



r a ia2 a 3 a 4 ( r i, r 2 , r 3 , r 4 ) -> W 2 lCt2 (n - r 2 ) (r 3 - r 4 ) + 

"ia 3 ( r i " r s) Wl 2ai (r 2 - r 4 ) + W 2 aiaA (n - r 4 ) W 2 ^ (r 2 - r 3 ) 



(96) 



In the same order of the loop expansion, the three point function will show up. The corre- 
sponding terms in kinetic part must cancel among themselves, as the potential term does 
not contribute. The local kinetic term yields the loop integrals of d 2 W 3 , whereas the non- 



local one yields VW 2 V'W 2 , integrated over d 3 r' with the Greens's function 

*9- 



(r-r') 
4n\r— r'\ 3 



The 





Fig. 7 

equation has the structure 

Now it is clear, that the solution of this equation for W 3 would be the same three point 
correlator, which one could obtain (much easier!) by direct iterations of the Navier-Stokes 
equation. 



= 
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The purpose of this painful exercise was not to give one more method of developing 
the expansion in powers of the random force. We rather verified that the loop equations 
are capable of producing the same results, as the ordinary chain of the equations for the 
correlation functions. 

In above arguments, it was important, that the loop functional belonged to the class 
of the abelian Stokes functionals. Let us check that our tensor area Anzatz 



C 



r Q dr 



belongs to the same class. Taking the square we find 



drp f dr' r a r' a 
c Jc ^ 



1 



2 f.. (lri y;."' < 



c 



dr'Jr — r') 2 



c 



(97) 



(9* 



where the last transformation follows from the fact, that only the cross term in (r — r') 2 
yields nonzero after double loop integration. 

Any expansion in terms of the square of the tensor area reduces, therefore to the 
superposition of multiple loop integral of the product of (r^ — rj) 2 , which is an example of 
the abelian Stokes functional. In the limit of large area, this could reduce to the fractional 
power. An example could be, say 





( 




\&[C] = exp 


B 


H 









/ 



\ 



(99) 



One could explicitly verify all the properties of the abelian Stokes functional. This example 
is not realistic, though, as it does not have the odd terms of expansion. In the real world 
such terms are present at the viscous scales. According to our solution, this asymmetry 
disappears in inertial range of loops (which does not apply to velocity correlators at 
inertial range, as those correspond to shrunk loops). 



B Matrix Model 



The Navier-Stokes equation represents a very special case of nonlinear PDE. There is a 
well known galilean invariance 

n Q (r, t) — > v Q (r — ut, t) + u a (100) 

which relates the magnitude of velocity field with the scales of time and space. f\ Let us 
make this relation more explicit. 

5 At the same time it tells us that the constant part of velocity if frame dependent, so that it better 
be eliminated, if we would like to have a smooth limit at large times. Most of notorious large scale 
divergencies in turbulence are due to this unphysical constant part. 
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First, let us introduce the vorticity field 

£*V = dpV,, - d u v^ (101) 
and rewrite the Navier-Stokes equation as follows 

v 2 

v a = vdpup a - vpujp a - d a w; w = p + — (102) 



This w is the well known enthalpy density, to be found from the incompressibility 
condition divv = 0, i.e. 

d 2 w = d a vptop a (103) 
As a next step, let us introduce "covariant derivative" operator 

D a = vd a - \ a (104) 

and observe that 

2 [D a Dp] = vujp a (105) 
2Dp [D a Dp] + h.c. = vdpujpu - vpujp a (106) 
where h.c. stands for hermitean conjugate. 

These identities allow us to write down the following dynamical equation for the 
covariant derivative operator 

D a = Dp [DaDp] - D a W + h.c. (107) 

As for the incompressibility condition, it can be written as follows 

D a D a ] = (108) 
The enthalpy operator W = ^ is to be determined from this condition , or, equivalently 

[D a [D a W}\ = [D a , D p [D a Dp]] (109) 

We see, that the viscosity disappeared from these equations. This paradox is resolved 
by extra degeneracy of this dynamics: the antihermitean part of the D operator is con- 
served. Its value at initial time is proportional to viscosity. 

The operator equations are invariant with respect to the time independent unitary 
transformations 

D a ^>S*D a S; S^S= 1 (110) 
and, in addition, to the time dependent unitary transformations with 

S(t) = exp (1-tuf, (Dp - Dj)) (111) 
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corresponding to the galilean transformations. 
One could view the operator D a as the matrix 



(i\ D <x\j) = J d 3 r^*(r>3 Q ^-(r) - 1 ^' k i {r)v a {r)'4) j {r) 



(112) 



where the functions ipj{r) are the Fourier of Tchebyshev functions depending upon the 
geometry of the problem. 

The finite mode approximation would correspond to truncation of this infinite size 
matrix to finite size N. This is not quite the same as leaving N terms in the mode 
expansion of velocity field. The number of independent parameters here is 0(N 2 ) rather 
then O(N). It is not clear whether the unitary symmetry is worth paying such a high 
price in numerical simulations! 

The matrix model of Navier-Stokes equation has some theoretical beauty and raises 
hopes of simple asymptotic probability distribution. The ensemble of random hermitean 
matrices was recently applied to the problem of Quantum Gravity ||, which led to a 
genuine breakthrough in the field. 

Unfortunately, the model of several coupled random matrices, which is the case here, 
is much more complicated then the one matrix model studied in Quantum Gravity. The 
dynamics of the eigenvalues is coupled to the dynamics of the "angular" variables, i.e. 
the unitary matrices S in above relations. We could not directly apply the technique of 
orthogonal polynomials, which was so successful in the one matrix problem. 

Another technique, which proved to be successful in QCD and Quantum Gravity is 
the loop equations. This method, which we are discussing at length in this paper, works 
in field theory problems with hidden geometric meaning. The turbulence proves to be an 
ideal case, much simpler then QCD or Quantum Gravity. 

C The Reduced Dynamics 

Let us now try to reproduce the dynamics of the loop field by a simpler Anzatz 



The difference with original definition (|5|) is that our new function P a (9) depends directly 
on 9 rather then through the function v a (r) taken at r a = C a (9). This is the d — *■ 1 
dimensional reduction we mentioned before. From the point of view of the loop functional 
there is no need to deal with field v(r) , one could take a shortcut. 

Clearly, the reduced dynamics must be fitted to the Navier-Stokes dynamics of original 
field. With the loop calculus, developed above, we have all the necessary tools to build 
this reduced dynamics. 




(113) 
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Let us assume some unknown dynamics for the P field 



Pa(d) =F a (9,[P)) (114) 

and compare the time derivatives of original and reduced Anzatz. We find in (|113|) instead 
of © 

l -<fdC a (6)F a (6,[P]) (115) 

Now we observe, that P' could be replaced by the functional derivative, acting on the 
exponential in (|1 13 ) as follows 



SC a (9) 



-ivP'(9) (116) 



This means, that one could take the operators of the Section 2, expressing velocity and 
vorticity in terms of the spike operator, and replace the functional derivative as above. 
This yields the following formula for the spike derivative 

D a (9, e) = -tuj e dcj) [ 1 - 1 Pl \ P'M = -iv jf_ i dix sgn(/i)P Q (9 + e(l + fi)) 

(117) 

This is the weighted discontinuity of the function P(9), which in the naive limit e — > 
would become the true discontinuity. However, the function P(9) has in general the 
stronger singularities, then discontinuity, so that this limit cannot be taken yet. 

Anyway, we arrive at the dynamical equation for the P field 

P a = vDfilfr - Vffrpa (118) 

where the operators V, D, Q of the Section 2 should be regarded as the ordinary numbers, 



with definition ( |1 1 7| ) of D in terms of P. 

All the functional derivatives are gone! We needed them only to prove equivalence of 
reduced dynamics to the Navier-Stokes dynamics. 

The function P a {9) would become complex now, as the right side of the reduced 
dynamical equation is complex for real P a (9). 

Let us discuss this puzzling issue in more detail. The origin of imaginary units was 
the factor of i in exponential of the definition of the loop field. We had to insert this 
factor to make the loop field decreasing at large loops as a result of oscillations of the 
phase factors. Later this factor propagated to the definition of the P field. 

Our spike derivative D is purely imaginary for real P, and so is our Q operator. This 
makes the velocity operator V real. Therefore the DQ term in the reduced equation ( 118|) 
is real for real P whereas the VQ term is purely imaginary. 

This does not contradict the moments equations, as we saw before. The terms with 
even/odd number of velocity fields in the loop functional are real/imaginary, but the 
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moments are real, as they should be. The complex dynamics of P simply doubles the 
number of independent variables. 

There is one serious problem, though. Inverting the spike operator D a we implicitly 
assumed, that it was antihermitean, and could be regularized by adding infinitesimal 
negative constant to D 2 a in denominator. This, indeed, works perturbatively, in each term 
of expansion in time, or that in size of the loop, as we checked. However, beyond this 
expansion there would be a problem of singularities, which arise when D 2 a {9) vanishes at 
some 9. 

In general, this would occur for complex 9, when the imaginary and real part of D 2 X {9) 
simultaneously vanish. One could introduce the complex variable 

e* = z; e-* e = <f dO = <f — (119) 

where the contour of z integration encircles the origin around the unit circle. Later, in 
course of time evolution, these contours must be deformed, to avoid complex roots of 
D 2 Q (9). 



D Initial Data 

Let us study the relation between the initial data for the original and reduced dynamics. 
Let us assume, that initial field is distributed according to some translation invariant 
probability distribution, so that initial value of the loop field does not depend on the 
constant part of C{9). 

One can expand translation invariant loop field in functional Fourier transform 

*[C] = J DQ5 3 (j> #<?(</>)) W[Q\ exp (i j> d9C a (9)Q a (9)^j (120) 
which can be inverted as follows 

5 3 (J d4Q{<f>)^ W[Q] = J DC*[C]exp (-i j> d9C a (0)Q a (9)^ (121) 

Let us take a closer look at these formal transformations. The functional measure for 
these integrations is defined according to the scalar product 

r rl9 

(A,B)=f-A(9)B(9) (122) 
which diagonalizes in the Fourier representation 

+oo 

A(9) = J2A n e md ; A_ n = A* n (123) 



+oo 



(A, B) = J2 AnB- n = A B + J2 <K + a'X; a' n = y/mA*, < = V2SA n (124) 
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The corresponding measure is given by an infinite product of the Euclidean measures 
for the imaginary and real parts of each Fourier component 

oo 

DQ = d 3 Q l[d 3 q> n d 3 ti (125) 
1 

The orthogonality of Fourier transformation could now be explicitly checked, as 

J DCexp (i J d9C a (9) (A a (9) - B a {9))) (126) 

= / d 3 C d 3 c' n dX exp (2m (c (A - B ) + £ c' n « - b' n ) + c» « - ) 

00 

= ^{A -B )\{5 3 K-b' n )5 3 K-b:) 
1 

Let us now check the parametric invariance 

9 -> /(2tt) - /(0) = 2tt; > (127) 

The functions C(9) and P(9) have zero dimension in a sense, that only their argument 
transforms 

C(0)->C (/(*)); P(0)->P(/(0)) (128) 
The functions <3(#) and P'(0) in above transformation have dimension one 

P'(0) - f'(9)P' (f(9)) ; Q(0) - f{9)Q (f(9)) (129) 

so that the constraint on Q remains invariant 

jd9Q{9) = jdf{9)Q{f{9)) (130) 

The invariance of the measure is easy to check for infinitesimal reparametrization 

f(9) = 9 + e(9); e(2n) = e(0) (131) 
which changes C and (C, C) as follows 

5C{9) = e{9)C'{9)- 5(C, C) = f -e(9)2C a (9)C' a (9) = - j -e'(9)C 2 a (9) (132) 
The corresponding Jacobian reduces to 

I- j>d9e'{9) = l (133) 

in virtue of periodicity. 

This proves the parametric invariance of the functional Fourier transformations. Using 
these transformations we could find the probability distribution for the initial data of 

P a {9) = -v f d4>Q a {4>) (134) 
Jo 
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The simplest but still meaningful distribution of initial velocity field is the Gaussian 
one, with energy concentrated in the macroscopic motions. The corresponding loop field 
reads 

9 [C\=exp(--fdC a (e) j>dC a {ff)f{C{6)-C{&))} (135) 
where f{r — r') is the velocity correlation function 

(v a (r)vp(r')) = (S^ - d a d p d- 2 ) f(r - r') (136) 

The potential part drops out in the closed loop integral. 

The correlation function varies at macroscopic scale, which means that we could ex- 
pand it in Taylor series 

f(r-r')^f -f 1 (r-r') 2 + ... (137) 
The first term f is proportional to initial energy density, 

\ (vl) = ^fo (138) 
and the second one is proportional to initial energy dissipation rate 

£ = -v (v a d 2 p v a ) = 2d(d - l)uh (139) 
where d = 3 is dimension of space. 



The constant term in (|137|) as well as r 2 + r' 2 terms drop from the closed loop integral, 



so we are left with the cross term rr' 

* [c] - ex P (-h jdc a (6) jdc a {e')c p {e)c^e')^ (uo) 

This is almost Gaussian distribution: it reduces to Gaussian one by extra integration 

%[C] -> const J d 3 o;exp (-a; 2 ^) exp Uijhu^ j dC^6)C v (6)^ (141) 

The integration here goes over all d ^ d ~ 1 ) = 3 independent a < (3 components of the 
antisymmetric tensor u a p. Note, that this is ordinary integration, not the functional one. 
The physical meaning of this cu is the random constant vorticity at initial moment. 

At fixed uj the Gaussian functional integration over C 

J dc exp (ijde (^Cp(e)p'p(6) + 2^ l u aP c'Mc P (e))) (142) 

can be performed explicitly, it reduces to solution of the saddle point equation 

I%B) = AvJhu^C'M (143) 
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which is trivial for constant u 

c a {e) = ^=u-^(e) (144) 

provided the matrix uj is invertible, which is true in general.^ We find the Gaussian 
probability distribution for P with the correlator 



(P a (9)P l3 (9')) = 2^ v // 1 ^sign(^ - 9) (145) 

Note, that antisymmetry of u compensates that of the sign function, so that this 
correlation function is symmetric, as it should be. However, it is antihermitean, which 
corresponds to purely imaginary eigenvalues. The corresponding realization of the P 
functions is complex! 

Let us study this phenomenon for the Fourier components. Differentiating the last 
equation with respect to 9 and Fourier transforming we find 



Au I — 

Pa,n,P/3,m) = ^-nmV fl^aB (146) 

m v 



This cannot be realized at complex conjugate Fourier components P a - n = P* n but 
we could take P a ,n = P a ~n and P a ^ n as real random variables, with correlation function 

Pa,nP/3,m) = — ^nmJYl^ap] U > (147) 

1 m v 
The trivial realization is 

4i/ / — 

Pa,n = —y flVapPp^ (148) 

with Pp n being Gaussian random numbers with unit dispersion. 

As for the constant part P a>0 of P a (9) , it is not defined, but it drops from equations 
in virtue of translational invariance. 



E W-functional 

The difficulties of turbulence are hidden in the loop equation, but they show up, if you 
try to solve it numerically. The main problem is that one cannot get rid of the cutoffs 
e, 5 — ► in the definitions of the spike derivatives. These cutoffs are designed to pick up 
the singular contributions in the angular integrals, but with finite number of modes, such 
as Fourier harmonics there would be no singularities. We did not find any way to truncate 
degrees of freedom in the P equation, without violating the parametric invariance. It very 
well may be, that this invariance would be restored in the limit of large number of modes, 
but it looks that there are too much ambiguity in the finite mode approximation. 



3 The integration over u> should be shifted towards complex plane to avoid such degeneracy. 
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After some attempts, we found the simpler version of the loop functional, which can 
be studied analytically in the turbulent region. This is the generating functional for the 
scalar products P a (9i)P a {9 2 ) 



W[S] = ^exp (-<fd0 1 <fd0 2 S(6 1 ,6 2 )P Q (0 1 )P a (e 2 ] 



(149) 



where, as before, the averaging goes over initial data for the P field. 
The time derivative of this W-functional 

w = -2(ldd 1 i de 2 s(e u e 2 )p a (9 1 )P a (e 2 )exp (- ld9 x <fde 2 s(e 1 ,e 2 )p a (e 1 )p a (e 2 



can be expressed in terms of functional derivatives of W by replacing 

6 



Pa(4>l)Pc 



(150) 
(151) 



for every scalar product of P fields, which arise after expansion of the spike derivatives 
(|TT71), (|I§), (|Sg) in the scalar product 



p a (ei)p a (o 2 ) = vPa{0x)Dp{e 2 )n Pa {e 2 ) - p a {9i)Vp{9 2 )n Poi {9 2 ) 

This equation has the structure 

w= <(d 2 es(e 1 ,e 2 ) \a 2 



(152) 



" 5 ' 




' 5 ' 




_5S_ 


W + A 3 


5S_ 





(153) 

where [X] stands for the k— degree homogenous functional of the function X(9\, 9 2 ). 

The operator D -2 is also the homogeneous functional of the negative degree k = —1. 
It can be written as follows 



with 



D- 2 (9,e)W\S]= drW[S + rU] 
Jo 

U(9x,9 2 ) = e~ 2 sgn(6> + e-0 x ) sgn(6> + e - 9 2 ) 



(154) 
(155) 



F Possible Numerical Implementation 

The above general scheme is fairly abstract and complicated. Could it lead to any practical 
computation method? This would depend upon the success of the discrete approximations 
of the singular equations of reduced dynamics. 

The most obvious approximation would be the truncation of Fourier expansion at 
some large number N. With Fourier components decreasing only as powers of n this 
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approximation is doubtful. In addition, such truncation violates the parametric invariance 
which looks dangerous. 

It seems safer to approximate P(9) by a sum of step functions, so that it is piecewise 
constant. The parametric transformations vary the lengths of intervals of constant P(9), 
but leave invariant these constant values. The corresponding representation reads 

N 

p a (o) = E M + 1) - p«(0) © (0 - ft) ; p{n + 1) = p(i), p(o) = o (ise) 

z=o 

It is implied that 6 = < 9i < 9 2 ■ ■ ■ < N < 2n. By construction, the function P{9) 
takes value p{l) at the interval Oi-i < 9 < 9\. 

We could take P(9) at the middle of this interval as approximation to p(l). 

p(Z)«P(ft); ft = ^(ft-i + ft) (157) 
As for the time evolution of angles ft , one could differentiate (|156|) in time and find 

p*(9) = E (p«(* + !) - p«(0) e (0 - ft) - E M + 1) - p«(0) <^ - W (iss) 

1=0 1=0 

from which one could derive the following approximation 

(P M (Z + 1)-J? M (0) A 

The extra advantage of this approximation is its simplicity. All the integrals involved 
in the definition of the spike derivative ( |117| ) are trivial for the stepwise constant P(9). 
So, this approximation can be in principle implemented at the computer. This formidable 
task exceeds the scope of the present work, which we view as purely theoretical. 
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